Three-Dimensional Elastic Compatibility: Twinning in Martensites 
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We show how the St.Venant compatibihty relations for strain in three dimensions lead to twinning 
for the cubic to tetragonal transition in martensitic materials within a Ginzburg-Landau model in 
terms of the six components of the symmetric strain tensor. The compatibility constraints generate 
an anisotropic long-range interaction in the order parameter (deviatoric strain) components. In con- 
trast to two dimensions, the free energy is characterized by a "landscape" of competing metastable 
states. We find a variety of textures, which result from the elastic frustration due to the effects of 
compatibility. Our results are also applicable to structural phase transitions in improper ferroelastics 
such as ferroelectrics and magnetoelastics, where strain acts as a secondary order parameter. 
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Strain plays a crucial role in many structural phase 
transitions, either as a primary order parameter (OP), 
e.g., in technologically important martensites [Q and 
shape memory alloys or as a secondary OP, e.g., in 
ferroelectric and magnetoelastic Q materials. More- 
over, self-consistent coupling of strain to either charge 
(e.g., high-temperature superconductors &) or spin (e.g., 
colossal magnetoresistance manganites g) degrees of 
freedom determines key physical properties. However, 
the physical degrees of freedom in a continuum elastic 
medium (such as described by Ginzburg-Landau models) 
are contained in the displacement field, u{r), even though 
it is the strain fields, e^, which appear in the free energy. 
Instead of treating the strains as independent fields, one 
must assume that they correspond to a physical displace- 
ment field, i.e., they are derivatives of a single continuous 
function. This is achieved by requiring that they satisfy 
a set of nontrivial (Saint- Venant) compatibility relations 
concisely expressed by the equation V x (V x e) =00 (in 
the absence of voids, dislocations, etc.). Indeed, ignoring 
the geometrical compatibility constraint and minimizing 
the free energy directly would lead to the erroneous result 
that non-OP strain components are identically zero, and 
that the OP responds trivially to perturbations, such as 
stress or local disorder, which is incorrect. 

We show for the first time that a consistent Ginzburg- 
Landau (GL) understanding of martensitic textures in 
terms of the OP-strain alone, and involving long-range 
OP strain-strain forces, is possible in three dimensions 
(3D). A 3D free energy expressed entirely in the order 
parameter strain variable (s), rather than the displace- 
ment field, provides a unified understanding of marten- 
sitic textures. We use 3D compatibility equations, link- 
ing the strain tensor components in the bulk and at in- 
terfaces, that induce anisotropic order-parameter strain 
interactions (masked in a conventional displacement-field 
treatment) . These two long-range bulk and interface po- 
tentials, together with local compositional fluctuations, 
drive the formation of global elastic textures. In contrast 
to 2D, the rich microstructure is the result of frustra- 



tion and competing metastable states in a complicated 
energy landscape. Here we will focus only on twinning 
for illustrative purposes. Other complex textures follow 
from our formalism and will be reported elsewhere. Our 
work is an essential step towards understanding results of 
High Resolution Electron Microscopy (HREM) and neu- 
tron scattering studies in alloys such as FePd and NiAl. 
In FexPdi_x, (Oil) twins preferentially grow at the ex- 
pense of (101) twins as the temperature is gradually low- 
ered to the martensitic transition temperature Q. In 
NixAli-x, martensite forms with a long period stacking 
structure involving (110) planes Due to the direc- 
tional aspects of the long-range interactions encoded by 
the underlying symmetry of the lattice, our analysis pre- 
dicts such twinning planes and directions. 

Model. As a specific example, we examine a cubic to 
tetragonal transition. In this case the symmetry-adapted 
strains can be defined in terms of the Lagrangian 
strain tensor components as the dilatation ei = 
(l/\/3)(eii -I- £22 + £33), the deviatoric (i.e., OP) strains 
62 = (1/V2)(eii -£22), 63 = (l/\/6)(eii-l-e22-2e33) and 
the three shear strains 64 = 2e23, 65 = 2ei3, eg — 2ei2- 
The elastic energy in 3D is given by ||l^ 

F = ^^(62,63) -I- -Fgra<i(Ve2, Ve3) + Fcs(ei, 64, 65,65), 

where the Landau free energy (Fig. 1), which provides 
the driving force for the transition, is 

Fl = Aiej + el) + Be^iel - iel) + C{e\ + e\f , 

and the Ginzburg (or gradient) contribution, which is re- 
sponsible for the interfacial energies, to the free energy 
is 
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The harmonic elastic energy contribution due to non-OP 
compression-shear (CS) strain components is 

Fcs^A^el + As{el + el + el). 

It is also possible to add a coupling, FcomposiVi '^Vi ^2, 63), 
of the OP strain to compositional fluctuations r/, which 
are responsible for tweed, but we do not consider such 
effects in the present paper. 
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FIG. 1. Countour plot of the Landau free energy in the 
martensite phase depicting the three degenerate energy min- 
ima corresponding to the three tetragonal variants with 
tetragonal axis along the x, y and z axis, respectively. Pa- 
rameters are A = -1.0, B = 0.83 and C = 0.04. 

The coefficients A = Ao{T — To) , B and C are related 
to the second, third and fourth order elastic constants, re- 
spectively, and can be obtained from experimental struc- 
tural data. The martensitic transition temperature is 
denoted To and the gradient coefficients g and h are de- 
termined from the phonon dispersion at the appropriate 
high-symmetry point in the Brillouin zone of the mate- 
rial. Ac and As denote the bulk compression and shear 
moduli, respectively. Here we do not consider the gra- 
dient energy contribution from the non-OP components 
ei, 64, 65, and eg, since their contribution is of secondary 
nature compared to the terms reported above. 

3D compatibility and analysis. We express the non- 
OP strains in terms of the OP strains (62, 63) using the 
three compatibility constraints £12,12 = £11,22 + £22,11, 
£23,23 = £22,33 + £33,22, £31,31 = £33,11 + £11,33- Analogous 
to the 2D case the constrained minimization of the 
compression-shear energy leads to an anisotropic long- 
range interaction between the OP strains in the bulk: 

pbuik ^ f d^rec,{a)U^c,,{\r-r'\)e^,{r') (a,a' = 2,3), 



where the kernels U22, U23, and C/33 can be expressed 
explicitly in Fourier space with each of the three kernels 
becomes a ratio of algebraic combinations of fc^,, ky and 
k^. All these combinations have the same order in the 
length, k, of the wavevector k. Hence, the kernels do not 
have any scale dependence. That is, the kernels only de- 
pend on the angles (j) and 9, defined in terms of spherical 
coordinates. The minima of each kernel occur at certain 
9 and (p values. These angles vary from kernel to ker- 
nel and the system, in general, experiences frustration. 
Figure 2 shows contours of the kernel U22 and U33 as 
functions of the two angles 9 and 0. The minima define 
a direction in 3D in which the compression-shear energy 
arising from the particular kernel is minimal. 

In contrast to 2D, Fes can be viewed as a weighted sum 
of the three kernels, where the fields 62 and 63 act as the 
weights. There is thus an interplay between Fl and the 
directional dependence of Fes to accommodate the frus- 
tration and choose minima, which in general do not have 
zero energy. Thus, there exist metastable minima cor- 
responding to different microstructures. In 2D there is 
no such interplay, as the directional minimum is inde- 
pendent of the amplitude of the field: There is only one 
minimum corresponding to zero energy and hence there 
is no frustration. Moreover, these long-range kernels in 
3D fulfill the full cubic symmetry properties associated 
with the transformation. In particular, Fes is invariant 
with respect to two-fold, three-fold, and four-fold rota- 
tion symmetries of the cube. We do not present here 
the mathematical structure underpinning this, but it is 
reflected in the frustration seen in 3D that is lacking in 
2D. 

The bulk contribution comes from a solid with periodic 
boundaries while the the effects of interfaces/surfaces 
such as habit planes (i.e., austenite-martensite or parent- 
product boundaries) may be included through a free en- 
ergy contribution, F^"^^ , where F^s = F^'^^^ + F^f^f and 

II 

= I d^kT{k,,ky,k,)[e2{kf + e3(k)% 

In analogy with the 2D case the function r(fc) has the 
form 



_ I kx i /Cy I I fcy Zb /C^: I I fcjs lb kx I _ 

for the [110] family of habit planes. 

We remark that the physics of twinning is quite dif- 
ferent in 2D and 3D. In 2D both the gradient and the 
bulk compatibility terms vary as ~ 1 /r^ and thus do not 
lead to a length scale. However, the surface compatibil- 
ity term in 2D is essential for introducing a length scale. 
It varies as ^ l/r which then competes with the com- 
bined bulk compatibility and gradient terms to give rise 
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to a twinning width. In 3D, the gradient terms behave as 
^ and F^""' ~ Thus, there is aheady a length 

scale present and F^^^f merely renormalizes the length 
scale arising from the bulk compatibility and gradient 
terms. We have numerically verified the above scaling 
relations. On the basics of the above arguments we ex- 
pect the twin width to scale as the square-root of the 
length in 3D to hold on the basis of above arguments. 
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FIG. 2. Three dimensional compatibility kernels U22 (up- 
per) and f/33 (lower) as a function of the colatitude 6 and 
the azimuth angles. Dashed contours surround minima, 
and dashed dotted surround maxima, while solid contours 
show intermediate levels. The parameters are As — 1.2 and 
A, = 2.4. 



Texture simulations. The various martensitic tex- 
tures that realize minima of the energy are found 
from random initial conditions and relaxational dy- 
namics for the OP strains. That is, ea{r) — 

where the time t is 
2,3. 

Depending on the details of the initial configurations, 
various twinning textures reflecting the metastability 



_ <5-F(ea,ei(ea),e4(ea),e5(eQ),e6(ea)) 
Sec^ (r) 

scaled with a characteristic relaxation rate and a 



emerge as ia vanishes. The simulations are performed 
for a regular cube, and periodic boundary conditions are 
applied in all six directions to obtain fully relaxed tex- 
tures in 62 (r) and 63 (r). 








X 



FIG. 3. (Color) Three dimensional twins in the (110) plane 
obtained from the time-dependent Ginzburg-Landau simula- 
tions, with representative parameters as given in Figs. 1 and 
2 and, h — 1 and g = 1. The scale of the twins is deter- 
mined by the competition between the interfacial energy and 
the bulk compression-shear energy, whereas the sharpness of 
the domain walls is determined by the Landau energy. 

We illustrate the simulation results obtained with rep- 
resentative (scaled) parameters. The red/blue/green 
color represent (62,63) positive/negative/zero OP strain 
values, respectively (the actual values are approximately 
determined by the minima of Fl as illustrated in Fig.l). 

Figure 3 shows the simplest twinning texture ((110) 
twins) that is uniform in the z direction. The strain 62 
changes from negative (minimum I in Fig. 1) to pos- 
itive (minimum II in Fig. 1) (i.e. blue to red). The 
strain 63 is not shown, as it is essentially constant with a 
small modulation at the twin boundaries. This particu- 
lar martensitic texture is the trivial extension of the 2D 
results since the OP fields are uniform in the z directions. 
Notice also that the directional dependence is the same 
as observed in 2D studies ]l^ . 

A richer texture is shown in Fig. 4, where the strain 63 
changes from negative (minimum III in Fig. 1) to posi- 
tive (minimum I or II in Fig. 1) (i.e. blue to red). The 
strain 62 would either be zero (green) when 63 is negative 
(minimum III) or positive/negative when 63 is positive. 
The figure clearly shows that twins with strain 62 can 
be oriented in two possible ways. This refiects the inver- 
sion symmetry of the kernels with respect to x, y. We 
note that a similar alternating twin microstructure has 
been implied by neutron and x-ray scattering studies in 
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layered high-temperature superconductors 
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FIG. 4. (Color) Three dimensional twins with two different 
orientations, (110) and (110), in adjacent planes. Both OP 
fields are shown, es in the upper panel and 62 in the lower. 

To illustrate the of frustration that is inherent in the 
compatibility relations, we explain the orientation of the 
twins in Fig. 4. We construct an effective kernel for the 
red and blue twins, respectively, using the appropriate 
values of the strains 62 and 63 in the twins. Figure 5 
depicts the minima of these effective kernels in the (f>-9 
plane. The dashed line shows the minimum correspond- 
ing to the blue twins and the solid line the minimum for 
the red twins. We note that the minima correspond to 
different <j> angles (45 + S and 45 — 6). The shift 6 depends 
on the parameters of the Landau free energy. Ideally, the 
red and the blue twins would not run parallel; however, 
the system accommodates these two conflicting directions 
by choosing the average angle of 45 degrees. This aspect 
of competing metastable states is a novel feature associ- 
ated with the cubic symmetry kernels that leads to a rich 
landscape of microstructures. 



FIG. 5. Minima of effective kernels illustrating inherent 
strain (orientational) frustration in 3D. 

Conclusion. We have derived a complete symmetry 
based, fully 3D model that describes the cubic to tetrag- 
onal structural transitions observed in many martensitic 
and shape- memory alloys. We obtained analytically the 
compatibility-induced anisotropic long-range potential in 
the OP (deviatoric strain tensor components) in the bulk 
and at interfaces. Unlike the 2D case 0,[2), there IS 
a length scale in the system from bulk compatibility 
alone; the surface compatibility potential also introduces 
a length scale akin to 2D. Many twin orientations are 
possible in 3D as a result of elastic frustration and the 
"landscape" (probably "rugged") of metastable energy 
states. Our model can be generalized straightforwardly 
to improper ferroelastics transitions by including symme- 
try allowed polarization (magnetization, etc.) nonlinear 
terms, and couplings to strain. Other symmetries can be 
handled within our approach. For example, we can study 
the cubic to trigonal (rhombohedral) transition in lead 
orthovanade, and NiTi- and AuCd-based shape-memory 
alloys using the three shear strains 64, 65, and eg as OP 
and with ei, 62, and 63 expressed in terms of the shear 
strains via compatibility. Image reconstruction methods 
of systematic experimental data from HREM and neu- 
tron scattering will allow a more quantitative comparison 
of strain patterns with those obtained from our model. 
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